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Abstract. The Painleve analysis introduced by Weiss, Tabor and Carnevale 
(WTC) in 1983 for nonlinear partial differential equations (PDE's) is an ex- 
tension of the method initiated by Painleve and Gambier at the beginning of 
this century for the classification of algebraic nonlinear differential equations 
(ODE's) without movable critical points. In these lectures we explain the WTC 
method in its invariant version introduced by Conte in 1989 and its application 
to solitonic equations in order to find algorithmically their associated Backlund 
transformation. A lot of remarkable properties are shared by these so-called 
"integrable" equations but they are generically no more valid for equations 
modelising physical phenomema. Belonging to this second class, some equa- 
tions called "partially integrable" sometimes keep remnants of integrability. In 
that case, the singularity analysis may also be useful for building closed form 
analytic solutions, which necessarily agree with the singularity structure of the 
equations. We display the privileged role played by the Riccati equation and 
systems of Riccati equations which are linearisable, as well as the importance of 
the Weierstrass elliptic function, for building solitary waves or more elaborate 
solutions. 



The Painleve property, one century later, ed. R. Conte, CRM series in math 
ematical physics (Springer- Verlag, Berlin, 1998) 31 March 1998 
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Chapter 1 

Introduction 



During the past thirty years, the interest for nonlinear phenomena has been 
growing in different fields of modern physics, such as optics, fluid dynamics, con- 
densed matter, elementary particle physics, statistical mechanics, astrophysics. 
Although the manifestation of those phenomena varies according to the different 
fields, they present a common feature in their mathematical description. The 
link comes from their description by nonlinear evolution equations (i.e. PDEs) 
whose solutions represent the propagation of waves with a permanent profile. 
Moreover, the analytical methods for solving them are directly inspired by the 
works of the famous mathematicians L. Fuchs, H. Poincare and P. Painleve, as 
explained in Conte contribution, whose content is assumed known. 

The propagation of a bell-shaped solitary wave on water has been approx- 
imately e xpla ined by the mathematical physicists J. Boussinesq Jl5| and Lord 



Rayleighlf f 4| only thirty years after its experimental discovery by Scott Russell 
in 1844 [117]. The full explanation was later given in 1895 by Korteweg and de 



Vries (KdV |74| ) who derived the nonlinear dispersive equation 

u t + u xxx + 3{u 2 ) x = 0, (1.1) 

possessing the two-parameter (fc, r) exact solution 

k 2 1 

u sw (fc,6») = — sech 2 -, sech= , 6 = k£ + r, £ = x - ct, c = k 2 . (1.2) 

y J 2 2 cosh s s v ' 



The name soliton was introduced by Zabusky and Kruskal in 1965 [138 when 



they solved the initial value problem for KdV equation (1.1) and discovered 



solutions describing the elastic collision of several waves (1.2). 

In these lectures we shall restrict our study by the method of singularities 
to nonlinear evolution equations possessing two different levels of integrabil- 
ity : complete integrability or partial integrability, including some chaotic PDE's 
which possess explicit analytic solutions in very special circumstances. 

Complete integrability means that 
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• either the nonlinear partial differential equation can be related to a linear 
partial differential equation by an explicit transformation, 

• or the equation passes the Painleve test and possesses the Painleve prop- 
erty (PP) for PDEs, i.e. : firstly, on every noncharacteristic manifold its 
general solution has no movable critical singularities in the complex plane 
of an arbitrary function tp(x,t); secondly, the PDE possesses an auto- 
Backlund transformation or is related by a Backlund transformation to 
another PDE possessing the PP (PDEs passing the " weak Painleve" test 
(I] and related by a hodograph transformation J27J to another equation 
possessing the PP are outside the scope of these lectures), 

• or the equation possesses solitary waves , A— soliton solutions for arbitrary 
A, an infinite number of conservation laws, bi-Hamiltonian structures, 
infinite-dimensional Lie algebras, . . . 

• or the equation satisfies the Ablowitz-Ramani-Segur (ARS) || || conjec- 
ture on the relationship of all its reductions to ODE's without movable 
critical points. 

More explicit definitions concerning the properties of this first class of equations, 
as well as classical examples will be given in chapter [2] 

Partial integrability means that some above listed properties are not satisfied 
(in particular the Painleve test may be satisfied only with some constraints 
on the function ip, or may never be satisfied whatever be ip, and the ARS 
conjecture is no more valid) but the equation possesses explicit analytic solutions 
like for instance : degenerate solitary waves, A— shock solutions, A— soliton 
solutions with A bounded p8j , or retains some pieces of integrability 

like degenerate Backlund transformations, a finite number of conservation laws 
, [58]]. For equations belonging to this second class, methods for finding- 
particular solutions, which must agree with the singularity structure of the 
equation, will be developed in chapter ||. 

Ch apte r || contains the main subject of our lectures : it is devoted to the 
WTC |13q| method and its extensions in the invariant version introduced by 
Conte [ |29[ , for finding algorithmically the auto-Backlund transformation of in- 
tegrable nonlinear PDEs. 
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Chapter 2 

Integrable equations 



We present here a few classical examples of nonlinear partial differential 
equations either explicitly related to linear partial differential equations or char- 
acterised by the properties of complete integrability mentioned in the previous 
chapter. 

2.1 Integration by direct linearisation 

Some equations can be linearised by an explicit transformation : 
Example 1. The Burgers equation 

u t + (u x +u 2 ) x = (2.1) 

is linearised into the heat equation p^ ] 

u = {Logip) x , ((cp t + ip xx )/ip)x = 0. (2.2) 



Example 2. The generalised Eckhaus equation 1 20 , [26| [78| 

iu t + u xx +(/3 2 \uf + 20e^(\u\ 2 ) x )u = O, (0,7) eft, (2.3) 
is linearisable into the Schrodinger equation for (3 cos 7 7^ 



iv t +U xx = Q,U = J-rl ^ e -(V2)tan 7 Log V with ^ = |„|2 ( 2 .4) 

V 2pcos7 y/tp 

and |u| 2 = — (Log«p) x . (2.5) 

2p cos 7 

If 7 = tt/2, the Kundu [[76], [nj gauge transformation u = ve % ^ v transforms the 
more general higher order nonlinear Schrodinger equation (HNLS) 

iu t + u xx + S\u\ 2 u + ([3 2 \u\ 4 + 2i(3\u\ 2 ) x )u = 0, (/?, 6) G ft, (2.6) 



G 



into the nonlinear Schrodinger equation (NLS) 

m + v xx + 5\v\ 2 v = 0. (2.7) 
The natural question is then : where do these miraculous transformations 



from u to another field come from? This will be answered in section 3.2.1 



2.2 Reduction to ODEs with the Painleve prop- 
erty 



Ablowitz, Ramani and Segur js], ||, Q, and McLeod and Olver [|90| 
tured a link between integrable NLPDEs and the Painleve ODEs 



107] 



conjec- 
for 



the integrable NLPDEs specially studied in these lectures, all known reductions 
to ODEs are singlevalued algebraic transfo rms of the Weierstrass or Painleve 
equations. Some of them are listed in Table 2.1. 



2.3 Construction of solitary wave solutions 

In the integrable case, the solitary waves sech and sech 2 are degenerate 
elliptic functions, obtained by imposing boundary conditions to the general 
solution of the ODE defining the travelling wave reduction. In the partially or 
nonintegrable case, the general solution of the reduction may not exist. One 
then looks for particular solutions, taking advantage of the singularity structure 
of the ODE by the method of subequations in chapter [|. 

Example 1. KdV 

The reduction u(x, t) = U(g), £ = x-ct,ol Eq. (Q yields the ODE 

(-cU + U" + 3U 2 )' = (2.8) 

After two integrations, this equation becomes 

-cU 2 /2 + U 3 + U' 2 /2 + K 1 U + K 2 = 0, (2.9) 

which identifies to the Weierstrass elliptic equation 

p' 2 = 4p 3 - 92p - .93, (92, 93) real constants, (2.10) 
u = c/6 - 2p(x -ct- x , c 2 /12 - K u K 2 /2 + K lC /12 - (c/6) 3 ). (2.11) 



The solitary wave (1.2) is found by imposing the boundary conditions {/(£) — > 0, 
U'(0 -> 0, E7"(f) -> 0, when |£| -> 00. Note that, for K x = K 2 = 0, equation 
( pTg| ) is a degenerate elliptic equation and 

p{x-ct-x ,c 2 /12,~(c/6) 3 ) = -(c/4)sech 2 (</c(x - ct- x )/2)+c/12. (2.12) 
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Example2. The generalised Tzitzeica equation 



u xt + ae u + a 1 e^ u + a e 2u = 0, fl^O 



(2.13) 



includes Liouville (ai 



«o 



0), sinh-Gordon (a x ^ 0, a = 0) or Tzitzeica 



[123, 124 1 (a 7^ 0, a! = 0) equations. It is polynomial in the variable v = e u 

w xt — v x Vt + av 3 + a\i) + do = (2-14) 

Its reduction (v, x, t) — > (V, £ = x — ct) can be integrated once 

- cV' 2 + 2aV 3 - 6KV 2 - 2a x V - a = 0, K arbitrary (2.15) 

and possesses the general two-parameter solution ]32^ 

aV = K + 2cp(£ - £ , (3X 2 + aai)/c 2 , (AK 3 + 2aa x K + a 2 a ) / (Ac 3 )) . (2.16) 

Moreover a linear superposition of two waves with opposite directions 

v(x,t) = Af(x - ct) + Bg(x + ct) (2.17) 

is compatible with the Tzitzeica equation by assuming that / and g satisfy the 
following second order ODE with constant coefficients 



f" = A x f + B x , g" =A 2 g 2 + B 2 

A particular solution of ( 2.13| ) for a x — is then 1 1 1 ] : 

ae u = 2c p(x - ct - x x ,g 2 ,K + a 2 a /(8c 3 )) 
-2c p(x - ct - x 2 , 52, K - a 2 a Q /(8c 3 )), 
c, xi, X2, 52, K arbitrary constants. 

Examples. NLS 



(2.18) 



(2.19) 



The reduction u(x, t) = p (^) e ^~ nt +^ of (|2.13l|) yields the coupled ODE's 

-cp' + 2<p'p' + ip"p = 0, 
p" + (n - (ip') 2 + ap')p + 2ep 3 = 

Equation ( [2.20 ) admits the integrating factor p 

tp' = c + K x /S, S = p 2 . 

Then Eq. ( [2.21 ) admits the integrating factor p' , hence 

S' 2 = -AeS 3 - AaS 2 + 8K 2 S - K 2 , a = n + c 2 /A, 

an elliptic equation for S with the general solution 



(2.20) 
(2.21) 

(2.22) 

(2.23) 



-a/(3e) - p(x -ct- x ,g 2 ,g 3 )/e 



(2.24) 



fj2 



8e(K 2 + a 2 /(6e)), g 3 = (2a/3) d + 8K 2 ae/3 + e l K{. (2.25) 
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The one-soliton solution is obtained for the values of K\ , K 2 making the Weier- 
strass elliptic function degenerate into a trigonometric function : 



92, 93) -> ai+a 2 sech fc£. 

This happens in two cases : 

1. ai = Ki = K 2 = 0, fc 2 = -a, p 2 = (k 2 /e) sech 2 fc£ 

2. aiKiK 2 ^0, ai = -(fc 2 +a)/(3e), p 2 = a x + (fc 2 /e) sech 2 fc£. 



(2.26) 



They respectively correspond for equation (2.131) to the three-parameter (c, fc, £0) 
solution ("bright" soliton) |l39| 



e>0: u = e- 1/2 ksech(k(x-ct-x ))e tcx / 2+ ^ k2 -^ 2 ^ t (2.27) 
and the four-parameter (c,k, K,xq) solution ("dark" soliton) ]140| 



e < : u 



{s)~ 1/2 [(fc/2) tanh(fc(a; - ct - x Q )/2) - i(K - c/2)] 



xc 



iA:x-2i[fc 2 /4+(K-c/2) 2 +K 2 /2]t 



(2.28) 



2.4 Conservation laws 

Definition. Given a PDE E(u;x,t) = 0, a conservation law is a relation 

T t + X x = 0, (2.29) 

where T and X, respectively called density and flux, depend on x,t,u and its 
derivatives. If the total variation of X in the interval a < x < b is zero, 
the quantity J T dx is a constant of the motion / called conserved quantity. 
"Integrable" PDEs possess an infinite number of conservation laws 128, 127. 
[|, For example, the first three conservation laws are : 



(a) for the KdV equation (yj) |136j |95j [72| 

Ti = u, Xi = 3m 2 + u xx ; 

T 2 = u 2 /2, X 2 = 2u 3 + uu xx - u 2 x /2- 

T 3 = 2m 3 — u 2 , X 3 = 9u 4 + 6u 2 u xx — 12uu 2 — 2u x u 3x 



(b) for the MKdV equation (^9J) |9J 

Ti = u, Xi = 2u 3 + u xx ; 

T 2 = u 2 /2, X 2 = 3u 4 /2 + uu xx - u 2 x /2; 



(2.30) 
(2.31) 
(2.32) 



(2.33) 
(2.34) 



T 3 = u 4 /4 - ui/4, X 3 = u b + u A u xx - Zu 2 u\ - u x u 3x /2 + t£j4; (2.35) 
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(c) for the sG equation fl2.118|) pi, Im 



Ti = u 2 x /2, X x = cosw; (2.36) 
T 2 = u x /A-u 2 xx , X 2 =ul cos it; (2.37) 
T 3 = 3u° - Ylu\u\ x + 16u 3 x u 3x + 12u\ x , X 3 = (2u x - 24^J cosu. (2.38) 



(d) for the NLS equation ( 2.131 ), we rep roduce th ree of the five conservation 
laws given by Zakharov and Shabat 13E, |l4Cf| for e = ±1 



h 
h 



\u\ 2 dx, I 2 



-oo 



(uu x — uu x ) dx, 



(K 



|m| 4 ) da: 



(2.39) 



OO 



/-t-oo 
(1 - \u\ 2 ) dx, I 2 = — (uu x - uu x ) dx, 

-oo J — oo 

(M 4 + K| 2 -l) dx (2.40) 



/a = 



(where u denotes the complex conjugate of u). 



For the Tzitzeica equation ( 2.62 ), Dodd and Bullough J43J first obtained two 
nontrivial conservation laws, then Mikhailov |94 gave a recursion formula for 
an infinite set of nontrivial polynomial conserved densities. 



2.5 Backlund transformations 



2.5.1 Definition 

A Backlund ||] transformation (BT) between two given PDEs 

Et(u; x, t) =0, E 2 (v;x',t') =0 (2.41) 

is a set of four relations (|^2| vol. Ill chap. XII) 

Fj(u,v,u x ,v x >,ut,v t >, ■ ■ ■ ;x,t,x',t') =0, j = 1,2 (2.42) 
x' = X(x,t,u,u x ,u t ,v), t' = T(x,t,u,u x ,u t ,v) (2.43) 

such that the elimination of u (resp. v) between (Fi,F 2 ) implies 

E 2 (v;x',t') = (resp. E\(u;x,t) = 0). 

In case the two PDEs are the same, the BT is called an auto-BT. 
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Backlund theory originates from the work of Lie and Backhand for the study 
of surfaces in differential geometry. The subject was subsequently developed by 
Goursat |57| and Clairin [ p4| . Backlund transformations re-present an extension 
of Lie contact transformations. They were first obtained for second order PDEs 
in two independent variables, linear in the highest derivatives (i.e. a special type 
of Monge- Ampere equation). 

For more details on BTs, the reader is advised to consult the book by Rogers 



and Shadwick 1 115 1 and the classical book of Goursat |5£ 



2.5.2 Examples : second order PDEs 

Burgers and heat equations 

Given the two equations 

E x = u t + (u x + u 2 ) x = 0, E 2 = v t + v xx = 0. (2.44) 
the two relations defining the BT are : 

Fi = v x —uv = 0, F 2 = v t + u 2 v + vu x = 0. (2-45) 

Indeed, the elimination of v (resp. u) yields the identities 

(F 2 /v) x - {Fx/v) t = Ei, v ^ 0, and F 2 + F hx + uF x = E 2 . (2.46) 

Liouville and d'Alembert 

Given the two equations 

E x = u xt - e" = 0, E 2 = v xt = 0, (2.47) 

the two relations 

F x = u x -v x + Xe (u+v)/2 = (2.48) 

F 2 ee u t + v t + {2/X)e^ u - v ^ 2 - 0, (2.49) 

where A is an arbitrary real constant called Backlund parameter, define a 
Backlund transformation as shown by the elimination of v (resp. u) 

Fx.t + F 2 , x - {l/\)e^- v ^ 2 Fx - {\/2)eS u+v ^ 2 F 2 = 2E t (2.50) 

Fx.t - F 2 , x + {l/\)e^- v ^ 2 Fx - (A/2) e ("+^/ 2 F 2 = -2E 2 . (2.51) 

Thus, the general solution of d'Alembert equation 

v = f(x)+g(t), (/, g) arbitrary functions, (2.52) 
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provides, by integration of the ODEs fl2.48D -( f^49| ), a solution of ( p.47| ) 



3 U = 2<p x <p t /<p 2 , = (A/2) / e f dx+ (1/A) / e^di 



which is the general solution. Travelling waves are built by the choice 



ip = cotanh(ax) — tanh(/3i) 



2a/?/ cosh 2 (ax - (3t). 



(2.53) 



(2.54) 



Sine-Gordon 

Given two solutions u and U of the sine-Gordon equation 

Ei = U x t - sin u — 0, E 2 = U xt - sin 17 = 0, 
the auto-Backlund transformation is defined by 
Fi = (u + f7) x -2Asin((u-{7)/2)=0 



(2.55) 



(2.56) 



F 2 



(it - J7) t - (2/A) sin((u + U)/2) = 0, A arbitrary constant, (2.57) 



as can easily be checked quite similarly to the Liouville and d'Alembert case, 
by elimination of U (resp.w) between these two relations 

Fi, t + F 2 , x + (1/A) cos((w + U)/2)Fi + A cos((u - U)/2)F 2 = 2E X (2.58) 
Fi, t - F 2 , x - (1/A) cos((u + U)/2)F X + A cos((u - U)/2)F 2 = 2E 2 . (2.59) 

Lamb [g0| built from ( 2.56] )-(2.57) infinite families of solutions, e.g. the N— soliton 
solution : at the first iteration, one starts from the solution U — ("vacuum"), 
and the integration of the ODEs ( 2.56)— (2.57) yields 



tan(u/4) = e Xx+x lt+s , S arbitrary constant, 
i.e. the one-soliton solution 

u x = 2A sech(A.T + X^t + 6), u t = 2A _1 scch(Ax + A _1 i + S). 

Tzitzeica 

For the Tzitzeica equation (Tzi) 

u —2u 

Uxt = e - e 



(2.60) 



(2.61) 



(2.62) 



there exists a complicated auto-BT [118, and another, much simpler one 
will be published soon |36|. We only report here the classical, well established 
results. 



12 



The Lax pair given by Tzitzeica |123| , 124 1 and rediscovered by Mikhailov 
33l B4| consists in the following matricial system of linear PDE's : 



(2.63) 
(2.64) 







[ d x <p 


■1 


V d t ^ J 








\ d x ip 


- 1 


V d t <p ) 








1 















e u 















1 





while the Moutard |97| transformation between two solutions U, u writes [ 124 1 : 

e u = -e u + 2 VxVt /ip 2 . (2.65) 

Since 1973, BTs have been found for PDEs of order greater than two. Dif- 
ferent approaches have been used for deriving those transformations : 



1. the method of Clairin B4L 82 



2. the method of differential forms developed by Wahlquist and Estabrook 



2 13fl and 



129j |l30j |46|), 

3. the method of bilinear transformations of Hirota |30| [32], [jjj , 

4. the method of gauge transformations developed by Boiti et al. 
Levi et al. M, 



In the last two methods, the BT results from the elimination of the wave 
function between the Lax pair and the DT. In next sections, these two main 
concepts of complete integrability are briefly recalled; then the principle of the 
method of gauge transformations is presented. Lax pairs and DTs are explicitly 
given for the PDE's of the AKNS scheme (KdV, M KdV, sine- Gordo n, NLS) 
and for some fifth order PDEs, respectively in sections 2.6.3 and 2.6.4 . But, let 
us first give some definitions. 



2.6 Darboux transformation and Lax pair 



2.6.1 Definitions 



Crum-Darboux transformation 

This transformation is a key in the theory of nonlinear integrable evolution 
equations for building soliton solutions and understanding their "asymptotically 
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linear" superposition rules. It is based on a result obtained by the French math- 
ematician Gaston Darboux in the special case of the Sturm-Liouville equation 
(also called Schrodinger equation in quantum mechanics). We briefly recall this 
old theorem and its generalisation due to Crum j40| . 

Theorem 1. (Darboux) The linear Schrodinger equation 

ip xx + {u + \)iP = (2.66) 

is invariant under 

V -> ^ = (9 X -%)V (2.67) 

u -> u = u + 2(Logip )xx (2.68) 

where ipo = ip[x, Xo) is an eigenfunction of ( [2.66] ) with parameter Ao- The 
essential point is that the new potential u depends only on -0o and not on 
tp. This transformation can then be iterated to obtain 

Theorem 2. (Crum) The function 

^ = w(^...^ N ) (2 - 69) 



where ijji, i/>2, ■ ■ ■ ipN are eigenfunctions of ( p.66| ) associated with parame- 
ters Ai, A2, . . . , Xn and th e sym bol W represents the Wronskian determi- 
nant, solves the equation ( 2.66| ) for the potential 



u = « + 2(LogW(^i,Va,..., M) xx - (2.70) 



Lax pair 



In 1968 Lax |83|| explained in a very transparent way the greater part of the 
result of Gardner et al. f35f by introducing the following operators 

L = -dl- u(x,t), A = -Adl - 6ud x - Su x (2.71) 

such that the KdV equation ([O]) may be represented in the following way 

d t L = [A,L] (2.72) 

called the Lax representation. Equation ( p. 72 ) expresses the compatibility be- 
tween the two partial differential equations of the system 

Llp = Xl/j 1 f l/j xx + (U + X)lp = (9 7-1) 

i^t^A^ J \ ^ t + (2u-4X)i> x -u x <tp = [ 
called Lax pair. This equivalence results from the identity 

4> xxt - i>txx = KdV(u)^ (2.74) 
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The system (2.73) is invariant under the Darboux transformation (2.67)-(2.68) 



with the compatibility condition 

(d t L)-4>=[A,L]i) (2.75) 
where (L, A) results from the substitution of u by u in (L, A). 

2.6.2 Backlund gauge transformation 

A general procedure to obtain BT for nonlinear PDE's derived as compat- 
ibility conditions between a given generalised Lax pair of operators was simul- 
tanuously considered by Boiti et al. and Levi et al. in 1982. It has provided 
new results for multidimensional nonlinear PDE's. Here, we only report the 
principle of the method. Let us consider the Lax pair 

4> x = Lip, Vt = Mtp (2.76) 

where tp is a N x N matrix as well as L, M which have a preassigned depen- 
dence on a matrix "potential" Q(x,t) and on a constant parameter A. The 



compatibility condition between the two equations of the system ( 2.76[ ) implies 
the following nonlinear equation 

L t - M x + [L, M] = 0. (2.77) 

To construct the BT for this nonlinear partial differential equation one has 



to consider two different systems of type (2.76) corresponding to two different 
"potentials", say Q(x,t) and Q(x,t) : 

i> x = L(Q(x,t);\)^ ^=M(Q(x,i);A)V (2.78) 
$ x = L(Q(x,t);X)$, 4> t =M(Q(x,t);X)$ (2.79) 

One assumes that the following generalised DT holds between the wave functions 
ip and tp : 

$ = Btjj (2.80) 

where B is a matrix function of Q,Q,x,t and A. The compatibility between 
( f2~8C| ) and the system (|2.78D-(f2~76|) gives the auto-BT 



B x = LB - BL, B t = MB - BM (2.81) 
By cross-differentiating these two relations one gets : 

(L t - M x + [L, M])B - B(L t - M x + [L, M]) = (2.82) 



which implies that if Q(x 7 t) satisfies the nonlinear PDE ( [2.77 ) then Q(x, t) 



satisfies the same equation. This exactly coincides with the definition of the BT 



previously given in section 2.5.1 
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Let us also mention the book of Matveev and Salle [B2| as a basis reference 
on Darboux transformation and its development in soliton theory. 

In the extension of Painleve analysis to NLPDEs |135[ | , if a PDE fulfills the 
necessary conditions of integrability ("Painleve test"), one tries to determine a 
Lax pair and a Darboux transformation relating two solutions of the same PDE 
in order to constructively prove the sufficiency of these conditions. A method 
(truncation procedure) leading to such a Lax pair and DT will be explained in 
section [3.2. 2| . In this formalism, the link with the notion of "general solution" is 
that the knowledge of the BT a priori allows to build wide classes of solutions. In 
one space dimension the "good" Lax pair of a given nonlinear PDE must depend 
on the solution of this equation and an arbitrary constant A. In the next section 
we show on examples how to derive the Lax pair and DT from the associated 
BT. In each case, it will be the aim of these lectures to show in chapter || how 
these two informations can be found algorithmically by singularity analysis. 

2.6.3 Examples : AKNS scheme 



Korteweg-de Vries 



Its conservative form is (1.1) and we define the potential form as 

u = w x , F(w) = w t + w XX x + 3w 2 x = 0. (2.83) 

Given two solutions w and W of ( 2.83| ), the auto-BT is defined by |82) 

(w + W) x = 2A - [w - W) 2 /2 (2.84) 
(w + W) t = -2(wl + w x W x + W2)-(w-W){w-W) xx , (2.85) 

where A is the Backhmd parameter. After changing variables w, W to W, Y = 
(w — W)/2, the gradient of Y is defined by the Riccati equations 

Y x = \-U -Y 2 , U = W X (2.86) 

Y t = (U x -(2U-4X)Y) x . (2.87) 

The transformation 

Y = d x LogiP (2.88) 

linearises these Riccati equations into one second order ODE and one first order 
PDE 

+ (U- = (2.89) 
V>t + (217 + 4A)V>* - (U x + G(t))tp = 0, G arbitrary function. (2.90) 

The Lax pair of KdV is defined by these two linear equations, which satisfy the 
compatibility condition 

Vw-V>t** = £TO, (2-91) 
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while the DT for KdV is defined by the a;- derivative of Eq. fl2.88| ) 

u-U = 2<9 2 LogV>. 

Modified Korteweg-de Vries 

Its conservative form is 

E(u) = u t + u xxx - 2a~ 2 (u 3 ) 
and we define the potential form as 

u = w x , F(w) =Wt + w xxx - 2a~ 2 wl 







Given two solutions w and W of (2.94), the auto-BT is given by 



(w + W) x = -2aAsinh((>- W)/a) 

(w + W) t = 8X 2 W x -4XW xx cosh((w-W)/a) 

+4a(2A 3 - XW 2 /a 2 ) sinh((w - W)/a) 

where A is the Backlund parameter. The change of variables 

(w, W) — (W, Y = e {w - w),a ) 

maps these equations into the two Riccati equations for Y, 

Y x = -2(U/a)Y + \(l-Y 2 ), U = W X 

Y t = 2A 1 r + s 1 (i + y 2 ) + c 1 (i-r 2 ) 

= (-4X17/ a + (2(U/a) 2 - 4A 2 + 2(U x /a))Y) x 



At 



U xx 
a 



U J 



4\ z 



u x u 2 , 

B 1 = -2A— , d = 2X— - 4A 3 . 
a a z 



o° a 

The compatibility condition of this "Riccati pseudopotential" Y is 

Y xt -Y tx = -(2/a)E(U)Y. 



(2.92) 



(2.93) 

(2.94) 

(2.95) 
(2.96) 

(2.97) 



(2.98) 
(2.99) 
(2.100) 

(2.101) 



(2.102) 



The Lax pair is obtained by linearising these two Riccati equations by the 
transformation 

Y = ipi/rfa (2.103) 



-U/a A \ />i 
A U/a)\i> 2 

B l - d -A l ) \ V^ 2 



(2.104) 
(2.105) 
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while the Darboux transformation is defined by 

u — U = ad x Log Y, 



which, by elimination of Y x with (2.98), is identical to 

u + U = aXiY- 1 - Y). 
The homographic transformation with a = U / a 

Y = X X /(l + a X ), 
maps the Riccati system ( ^.98 )-( |2.99| ) into the simpler form : 
Xx = 1 + (S/2) X 2 



Xt = - 
S = 2 



-C + C xX -(l/2)(CS + C xx ) X 2 

2 \ I 



a 



X 



c 



u x 



2X Z 



We shall see that the relation between the two functions S and C 

S - C + 6A 2 = 



(2.106) 



(2.107) 



(2.108) 



(2.109) 
(2.110) 

(2.111) 



(2.112) 



corresponds to the singular manifold (SM) equation of the KdV equation [135 
and can be found algorithmically [ 1 10 1 when one performs the Painleve analysis 
of the MKdV equation. In the variable 



f = a(Y-l)/(Y + l), 
the system (p.98|)-( |2.99D and the DT fl2.106p or ( |2.107| ) become 

u-U - 2a 2 f x /(a 2 - f), 

u + U = -Aa 2 Xf/(a 2 - f), 

a f x = -{U/a)(a 2 - f 2 )-2Xaf 

af t = A 1 (a 2 -f 2 ) + B 1 (a 2 + f 2 )-2C 1 af. 



(2.113) 



126, 12 



(2.114) 
(2.115) 
(2.116) 
(2.117) 



Sine-Gordon 



E(u) = u xt - smu = 0. (2.118) 
Given two solutions u and U of ( ^.118 ), the auto-BT is given by [p0[ 

(u + U) x = -4Asin((u-t7)/2) (2.119) 

{u-U) t = -A- 1 sin((u + U)/2), (2.120) 
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where A is the Backlvmd parameter. The change of variables 

(«, U) -> (U, Y = e -H»-v)/*) 

maps these equations into the two Riccati equations for Y 

Y x = iU x Y + A(l - Y 2 ) 

Y t = ((1 - Y 2 ) cost/ + + Y 2 ) sin J7)/(4A). 

The compatibility condition of the Riccati pseudopotential Y is 

Y xt - Y tx = iE(U)Y. 

The Lax pair is obtained by linearising the Riccati system 



i>2 



1p2 



iU x /2 A \ fipi 
A -iU w /2 J U 2 



1 

4A 



while the Darboux transformation for sG is defined by 

u-U = 2iLogF, 
(u + U) x = 2i\(Y- 1 -Y). 



(2.121) 

(2.122) 
(2.123) 

(2.124) 
(2.125) 

(2.126) 

(2.127) 
(2.128) 



The ho mograp h ic tran sformation ( 2.108 ) wit h a = —iU x /2 maps the Riccati 
system ( ^.122| )-( ^T2^ ) into the simpler form ( ^.109|) -( |2T10| ), with 



S = -iU xx + Ul/2 - 2A 2 , C = -e lU /(A\ 2 ). 
The relation between S and C 

S + C xx /C - (1/2){C X /C) 2 + 2A 2 = 



(2.129) 



(2.130) 



represents the SM equation obtained by Conte p9[ when performing the invari- 
ant Painleve analysis of the sine-Gordon equation. 

Nonlinear Schrodinger 



E(u) = iu t + u xx + 2e\u\ 2 u = 0, e = ±1. (2.131) 
Given two solutions u and u' of (|2.131 ), the auto-BT can be written as [^3[ 



73 B6 



(u + U) x = (u-U)y/4\ 2 -e\u + U\< 



(u + U) t = i(u-U) x ^AX 2 - e\u + U\ 2 

+ie{u + U)(\u + U\ 2 + \u - U\ 2 )/2. 



(2.132) 
(2.133) 
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The extension to NLS of the transformation ( 2.11 5| ) is 



u+U=-AXf/(l + e\f\ z ). (2.134) 



Therefore, the change of variables (it, U) — > (U, /) transforms (2.132) into 

-fx + ef 2 J x - (1 - e\f\ 2 )(U(l + e\f\ 2 ) + 2Xf). (2.135) 

The elimination of / between this equation and its c.c., assuming 1 — |/| 4 ^ 0, 
provides 

f x = ~2Xf-U-eUf 2 . (2.136) 

while the t— part is 

ft = (XU + U x ) + [eUU + A 2 )/ + {XU - U x )f 2 (2.137) 
with the identity 

fxt- ftx = E + eEf 2 . (2.138) 
Equations ( p.!32| ) and ( |2.134[ ) imply 



U -C/ = 2(/ K -/ 2 /J/(l-|/| 4 ). (2.139) 

In all the above examples (KdV, MKdV, sG, NLS), the DT is defined with 
one (for KdV) or two (for the others) entire functions ip. This distinction is the 
only relevant feature needed to obtain in an algorithmic way the Lax pair by 
methods linked to the singularity structure of these equations. 

2.6.4 Higher order KdV-type equations 

Among the fifth order nonlinear evolution equations 

U t + (u xxxx + (8a - 2(3)uu xx - 2(a + (3)u 2 x - (20 /3)a(3u 3 ) x = (2.140) 

only three cases are integrable : 

13 1 a = -1 : ut + (u xxxx + 10auu xx + 20a 2 u 3 /3) x = (2.141) 

(3/a = -6 : u t + (u xxxx + 20auu xx + 10cm 2 . + 40a 2 u 3 ) x = (2.142) 

0/a = -16 :u t + (u xxxx + 40auu xx + 30au 2 x + 320a 2 u 3 /3) x = (2.143) 

respectively named Sawada-Kotera (SK) or Caudrey-Dodd-Gibbon |l20| , pl| , 
Lax's 5th order KdV (KdV5) || and Kaup-Kupershmidt (KK) §§. 
Their respective Lax representation ( 2.72 ) is [^3[ ^3lJ 



(SK) a = 3, L = d 3 + 6ud x , (2.144) 
A = 9d x + 90udl + 90u x dl + (60u xx + I80u 2 )d x ; 
(KdV5) a = 1/2, L = d 2 + u, 

A = Wdl + 40udl + 60u x dl + {50u xx + 30u 2 )d x + 15u xxx + 30uu x ; 

(KK) a = 3/4, L = d 3 x + 6ud x + 3u Xl (2.145) 

A = 3{3dl + 30ud 3 + Vbu x d 2 x + (35u xx + 60u 2 )d x + 10u xxx + 30uu x ). 
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We discard the equation ( [2.142 ) for it has the same second order scattering 
problem as the KdV equation, and we restrict to the two equations Q2.141 )- 
( [2.143 ) possessing two different third order scattering problems Lip — \ip. 
The SK equation possesses the Darboux transformation 11£] 



u = U + d 2 x Log V> 
while for the KK equation this transformation writes J8lj 

u = U+ (1/2)% Log cp, ip = ^ xx - (1/2)^ + 3C/V 2 . 



(2.146) 



(2.147) 



In the notation w x = u,W x = U, the x-part of the BT for SK |4J, |119[ is 

(w - W) xx + 3(v - W)(w + W) x + {w- Wf = A, (2.148) 
while for the KK equation it writes 

(w ~ W) xx + 3(w - W)(w + W) x - (3/4)(w - W) 2 x /(w -W) + (w- W) 3 = X. 

(2.149) 

This last expression was obtained for the first time by Rogers and Carillo [ |116| 
in the particular case A = 0. 
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Table 2.1: Some reductions of a PDE to an ODE and their solutions. The PDEs 
E(u,x,t) = (K dV, M KdV, sG, Bq, NLS, T zi ) are respectively defined by the 
equations Q, ( ^93| ), ( |2.118[ ), ( p3(i| ), fl2.13l| ), ( ^62| ). The reduction to an ODE 
for U(£) is defined by the two expressions of u in terms of (U,x,t) and of £ in 
terms of (x,t). The letter K, with or without subscript, denotes an arbitrary 
constant. Last column indicates the elementary function (p, (P1)-(P6)) whose 
general solution of the ODE is a singlevalued algebraic transform. 



1 I 1 1 


u 




i \ 1 \ 1 ' 


P, (Bn) 


KdV 


U 


x — ct 


U' 2 + 2C/ 3 - cU 2 
+ 2KxU + 2K 2 = 


P 


KdV 


U - Xt 


x + 3Xt 2 


U" + 3U 2 
— A? + K = 


(PI) 


MKdV 


U 


x — ct 


U' 2 -U 4 - cU 2 
+ K 1 U + K 2 = 


P 


MKdV 




xiM)- 1 / 3 


U" - 2U 3 
-£U + K = 


(P2) 


sG 


—i Log U 


x — ct 


cU' 2 + U 3 
+ KU 2 + U = 


P 


sG 


—i Log U 


xt 


u" - u' 2 /u + u'fi 

+ (1 - U 2 )/{20 = 


(P3) 


Bq 


U 


x — ct 


(U"/3) + U 2 
+c 2 U + K 1 £ 
+ K 2 = 


P, (PI) 


Bq 


2(U' + ^-t 2 ) 


x-t 2 


U" 2 + 4U' 3 
+12{£U' - U)U' 
+ K 2 U' + K 3 = 


(P2) 
[113, H 


Bq 


(w-em/t 


xt- 1 / 2 


U" 2 /2 + U' 3 
-(9/8)(U-£U') 2 

+K 1 {U-iU') 
+ K 2 U' + K 3 = 


<g 


if 


NLS 


(2.20) 


x — ct 


(2.23) 


P 


NLS 


e i(xt-4t 3 /3)jj 


x-t 2 


U" + 2eU 3 
- 2£U = 




P2) 


NLS 




xt- 1 / 2 


4C/" 2 + 4eU' 3 + KU' 
+ - U) 2 /A = 




Tzi 


hogU 


x — ct 


- cU' 2 + 2aU 3 
+ KU 2 - a = 


P 


Tzi 


hogU 


xt 


{{U'/Uy + aU 
+a U- 2 = 
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Chapter 3 

Painleve analysis for PDEs 



The WTC extension [ 135 1 of Painleve analysis to partial differential equa- 
tions consists of two parts 

1. generation of necessary conditions (Painleve test) for the absence of mov- 
able critical singularities in the "general solution" , 

2. explicit proof of sufficiency by finding the transformation which linearises 
the PDE or yields an auto-BT or a BT to another PDE with the PP. 

The methods relative to both parts are different. 

In the first part, for every noncharacteristic manifold (v?(x, t) — 0,(p x 
0), one tests the existence of all possible local representations of the "general 
solution" by a Laurent series in the neighbourhood of tp = 0. This test may 

• pass whatever be ip; the PDE may then have the PP, 

• fail whatever be ip; this is typical of chaotic PDEs, 

• pass with some constraints on (p; then there exists particular Laurent series 
and the PDE is called "partially integrable" . 



In the second part, the Weiss truncation procedure [ 131 1 , using only the 
singular part of the Laurent series, may yield constructive results like 

• the linearising transformation or the BT, in case the PDE passes the Pain- 
leve test for every tp, 

• particular solutions, necessarily compatible with the singularity structure 
of the PDE, in case the Painleve test is conditionally or not satisfied (see 
chapter 0). 



23 



3.1 Necessary conditions (Painleve test) 



Contrary to the case of ODEs, the singularities in the complex domain of 
(x,t) are not isolated. Given a PDE E(u,x,t) = of order N polynomial in 
u and its partial derivatives (maybe after a preliminary change of variables), 
we consider the associated equation p(x, t) = of the movable SM and an 
expansion of u and E as a Laurent series in \ m the neighborhood of ip = 0. 
We distinguish between ip and the expansion variable x an d only require x to 
vanish as p 



u{x,t) = ^ Uj (x,t) X ^ p , E(u,x,t) = ^E 3 {x,t) X 3+ \ (3.1) 

j=0 j=0 

where (p, q) are two negative integers with q < p — 1, and (lij, Ej) the Laurent 
series coefficients. The result of the Painleve test (necessary conditions) is inde- 
pendent of the explicit expression for x but some particular choices are better 
than others during the second part (sufficient conditions) when one looks for 
the Lax pair or tries to linearise the equation. 

The main choices (gauges) for the expansion variable x are 



i) WTC gauge [ 135 1 x = <p, hence coefficients {uj,Ej) rational in the 



derivatives Dp of p 

• ii) dimensionless WTC gauge x — p/Px, hence coefficients (uj,Ej) ratio- 
nal in the derivatives Dp of p of homogeneity degree zero, 

• hi) Kruskal gauge J68[ x = x ~ /(*)) / arbitrary, hence coefficients (uj, Ej) 
independent of x and rational in the derivatives of /. This is the simplest 
choice for the test, but it cannot be used to obtain the Lax pair or partic- 
ular solutions. 



iv) Contc gauge |p9[ y = p/(p x - PxxP/{2px)) ~ v ~>o p/p x , hence coef- 
ficients (uj,Ej) rational in the derivatives of p invariant under the group 
of nomographic transformations 

ip — > (ap + b)/(cp + d), (a, 6, c, d) arbitrary complex constants. 

In this last case, the Riccati system satisfied by x is 

Xx = l + (S/2) X 2 (3.2) 

Xt = -C + C xX ~(l/2)(CS + C xx )x 2 (3.3) 

2((xr 1 )^-(x^ 1 )0 = S t + C xxx + 2C x S + CS x = (3.4) 

with 

S = {<p; x} = {p xx /ip x ) x - (l/2)(p xx /ip x ) 2 , C = -p t /ip x . (3.5) 
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The transformation x = ip/i>x linearises this Riccati system into 

?P XX + (5/2)V - (3.6) 
ipt + Cipx — {C x /2 + g(t))ip = 0, g arbitrary function. (3.7) 

This choice of gauge is equivalent to the expansion of (u, E) as 

+ oo +00 

u = 5> j (VM0 i+p , E = Y, E iW^) j+q , (3-8) 

j=0 j=0 

where the function ip satisfies a second order linear ODE in the x variable. 
To obtain the couples (uo,p) one substitutes in the polynomial PDE 

u - «oX p , Xx -» 1, Xt -» -C, £)u -» « ^(X P ). (3.9) 

One then determines the balance between the different terms of this polynomial 
expression. Each different solution (uo,p) defines a family. For every j > 1 the 
recurrence relation determining tij is 

Vj > 1 : PC7>j - Qj^Ufc.Dufc, fc G [0, j - 1]}), (3.10) 

where P is a polynomial of degree at most N. 
The main requirements of the Painleve test are 

• the zeros of P (Fuchs indices, also named Painleve resonances) are distinct 
integers, 

• for every index i and every ip, the compatibility condition Qi = holds. 

3.2 Methods for proving sufficiency 

One distinguishes two main methods : 

1. the singular part transformation which may provide the explicit transfor- 
mation linearising the nonlinear PDE. If this is not the case, the transfor- 
mation may yield an equation in a form more convenient than the original 
one to search for explicit solutions, 

2. the truncation procedure of Weiss and its extensions for obtaining the BT 
and thus proving that the nonlinear PDE possesses the PP. 

3.2.1 Singular part transformation 

The method consists of transforming the PDE for u into an equation for 
by the nonlinear transformation 
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u = V Log tp, 



(3.11) 



where T> is the singular part operator associated with one of the families defined 
in the Painleve test. 

Example 1 (linearisation). Burgers equation 

u t + u xx + (u 2 ) x = : u = (p. x (p~ l (3-12) 

u = PLog^ = d x Log if ](pt + (p xx +K(t)ip = 0, K(t) arbitrary function. (3.13) 
Example 2 (linearisation). Liouville equation 

Vxt - e v — 

e v = u, uu xt - u x u t - u 3 — : u — 2(tp x tp t ^ 2 - Pxt^ 1 ) (3.15) 



u = VLogip = -29^(Log<^), ip xt = 0. 
Example 3 (linearisation) . Eckhaus equation |2(], |2^, Q 

iu t + u xx + q r (\uf + 2a(\u\ 2 ) x )u = 0, a 2 = l/q r , q r £ TZ. 



(3.14) 
(3.15) 
(3.16) 



(3.17) 



In the variables (w, 6) defined by 
equivalent to the system 



axgu,w x 



2 , the equation (3.17) is 



- o— , t = - 2 f- ) - ~-L + q r ( w 2 x + 2aw xx ) 3.18 

2 w T 4 w x wi 4 u>;£ 



1 ^ l^xxx 

4 V ir, 



1 utf 



whose compatibility condition is 

x t ~ Qtx = (wttwl + w xx w 2 )/2 + {w xxxx w 2 x + wl x )/2 - w x w xx w xxx 

- w t w x w xt + 2q r (w x w xx + aw x w xxx ) = 0. (3.19) 

Under the transformation w — {a/2) Log</? defined by the singular part operator 
T> of the equation for w, these three equations become 

9 X = -<pt/(2<p x ), 8 t = Vxxx /(2 Vx ) - y 2 xx /(4v 2 x ) - VtK^l) (3-20) 
9 xt - 9 tx = (<pu<fl + ¥xx<Pt)/ 2 + {<Pxxxx<fil + <pU/2 

— fxfxx^Pxxx — WPxfxt = 0. (3-21) 



The t hree e quatio n s (|3.20| ) , ( |3.20| ) , ( p.2l[ ) are deduced from the three previous 
ones (3.18), (3.18), (3.19) by the following simple operation : change w to <p and 
assign q r to zero. Thus the transformation has linearised the Eckhaus equation 
( |3.17| ) into the Schrodinger equation 



wt + Vxx = 0, ip x = \v\ 2 , axgv = argw. 



(3.22) 
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Because of the conservation of the phase, one finally has 



u = (Va/2>/y J \is\ 2 dx, \u\ z = (a/2)d x Logtp. (3.23) 

Example 4 (bilinearisation). Korteweg-de Vries equation 

u t + u xxx + 3{u 2 ) x = : u = -2<p x ip~ 2 + 2ip xx (p~ 1 (3.24) 

u = Vhog V > = 2d 2 x hog V >- (D x D t + D x ){ip ■ ip) = Q (3.25) 

The transformed equation, quadratic in ip (see the numerous papers of Hirota 
[ [59"[ |6C[| for the definition of the bilinear operators D x , D t ) is convenient to look 
for iV— soliton solutions, auto-BTs, Miura transformations. 

3.2.2 Weiss method and its limitations 

If a nonlinear PDE passing the Painleve test is not linearisable, the idea of 



Weiss [135, 131] is that the principal part of this local Laurent series contains 
all the information for proving that the PDE possesses the Painleve property 
through the knowledge of its BT (i.e. its DT and Lax pair). This method consists 
of truncating the Laurent series for it and E(u) to their nonpositive powers in 

X 

UT=Y, u iX j+P , E T = Y1 E oX j+q (3.26) 

3=0 3=0 

and identifying to zero the coefficients Ej of the x~ polynomial \~ q Et(x)- 
Equations Ej = for j = 0, . . . , — p determine the p + 1 coefficients uj as equal 
to those of the infinite expansion. After replacement of Uj by these values, the 
remaining equations are 

^■(£V,O = j S {— p + 1, . . . , —q} (3.27) 
j =/= compatible indices , i = indices S {0, . . . , —p} 

In the Conte gauge, the coefficients Uj , Ej depend on the derivatives of ip through 
the homographic invariants (<S*, C) and their derivatives. As the variable x = 
ipx/ip satisfies a Riccati equation one can connect the monomial (ip x /tjj) n with 
the derivatives (Log^)^, (j < n 6 JV + ) and show that 

-p-i -p 

J2 Uj(S,C)^/^ x y + P=J2^(S,C)(Log^) jx + f(S,C). (3.28) 

3=0 3=1 



Then the relation 



Ur -fi = 5^fi i (S,C0(LogV)ix = I>LogV, (3.29) 

3=1 
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where V is the singular part operator and u = U- P (S,C) + f(S,C), defines 
a Darboux transformation if E(u) = 0. For this reason, we call equations 
( |3.27 ) Painleve- Darboux equations. The elimination of the arbitrary functions 
Ui among this set must produce only one "independent" equation 



F(S, C) = 



(3.30) 



called the singular manifold equation, modulo the ever present link between S 
and C given by equation (|3.4|). 

The next step consists of finding a parametric representation for equation 
( |3.3C| ) under the form (S, C) depending on a function U and an arbitrary con- 
stant A such that the cross-derivative condition (3.4) is identical to the original 
equation E(U) = for U. If this is indeed the case and if U can be identi- 
fied with u, the truncation will provide the DT (as a consequence of equation 
(p9|)), the Lax pair (as consequences of the linear system (|3.6|)-(|3~7|) and the 
parametrisation of S and C) and thus the BT. 

This method only succeeds for a few equations, like KdV |135|, KdV5, AKNS 



[ |98|, all belonging to the same hierarchy. Let us describe it for the KdV equation 
This equation, which passes the Painleve test, admits the single family 
u ~ — 2\~ 2 with Fuchs indices —1, 4, 6 



u = -2 X - 2 + (C - 45)/6 - (1/6)(C - S) xX + 0( x 2 ). 



(3.31) 



The algorithmic results of the Painleve analysis for KdV are given in Table 
1[ They yield the SM equation 

(3.32) 



C — S* + 6A = 0, A = arbitrary constant. 
Its parametric representation 

S = 2(U + X), C = 2(U-2X) 



provides the second order linear system (|3.6|) — (p.7|) 



^Jxx + (U + X)ip = 

ipt + 2(U - 2X)tf, x - U x ^j 







(3.33) 



(3.34) 
(3.35) 



satisfying the cross-derivative condition ip xx t — iptxx = 2 KdV(U)ijj = 0. The 
map between two solutions of KdV coming out of the truncation is 



u T = 2{Logi>) xx + {C + 2S)/6 = 2(Log^) ra + U 



(3.36) 



Thus, the Weiss truncation yields both the Lax pair (2.89)— ( 2.90) and the DT 
(|j|) of the KdV equa tion. The auto-BT Q2.84| )-( ^85j is obtained by su bsti- 
tution of the DT fl3.36| ), i.e. ip x /ip = (w - W)/2, into the couple ( ^.34[ ) and 
( 3.35 ) (notation iiy = u = w x , U = W x ). 
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It happens that, for other equations possessing either one family of mov- 
able singularities or several families with nonopposite residues like Boussinesq, 
Sawada-Kotera, Hirota-Satsuma ]63| equations, the parametrisation of (S, C) 
yields a condition ( |3.4D for U different from the original equation, this defines 
a transformation between ut and U , called Miura transformation, obtained by 
the elimination of (x, S, C) between the four equations : ut — the truncation, 
the two equations of the parametric representation (5, C) — f(U) and anyone 



of the two (nonindcpendent) equations (3.2), (p.3|). Then in o rder to obtain the 



auto-BT, one requires that the function tp in equation ( p. 29] ) satisfies a linear 
third order system whose coefficients are to be determined as functions of A and 
another solution U of the analy-sed PDE linked to ut through the Darboux 
transformation. 

3.2.3 Method for third order Lax pair 

Let us denote (a, b, c, d, e) the five unknown coefficients defining a third order 
linear system for tj; 

4>xxx = aip x + bip (3.37) 
ipt = cip xx + dtjj x + &ip (3.38) 

whose compatibility condition is 

{ipijxxx ~ (ipxxx)t = X ip + Xxip x + X 2 i> xx = (3.39) 

Xq = —bt — ae x + e xxx + b xx c 

+3bc xx + 3b x c x + 3bd x + b x d = (3.40) 

X\ = —at + 3e xx + 2b x c + a xx c + d xxx + 3ac xx + 2ad x (3.41) 

+3>a x c x + 3>bc x + a x d = (3.42) 

X 2 = (2ac + c xx + 3d x + 3e) x =0. (3.43) 

In t he tw o indep endent components Z\ = ip x /ip,Z2 = i^xx/ip, the linear 
system ( 3.37| )-( 3.38 ) is equivalent to the projective Riccati system 



Zx,x = {-Z 1 )Z 1 +Z 2 (3.44) 

Z 2 ,x = (-Z 1 )Z 2 +aZ 1 +b (3.45) 

Z ht = (-dZ 1 ~cZ 2 )Z 1 + (ac + d x )Z 1 + {c x +d)Z 2 +e x + bc (3.46) 
Z 2 ,t = (~dZx - cZ 2 )Z 2 + (2ac x + a x c + bc + d xx + ad + 2e x )Z 1 

+ {c xx + 2d x + ac)Z 2 + 2bc x + b x c + bd + e xx . (3.47) 

The determining equations for the coefficients (a, b, c, d, e) of the Lax pair 
are generated by the expansion of Et = E(ut) on the basis (Zi, Z 2 ) 

e t = y, Ci ^ z ™' ( 3 - 48 ) 
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Ci :m = Ci, m {a, b, c, d, e, U) = 0. 



(3.49) 



In case the solution of the determining equations does not lead to the ex- 
pected solution, for a reason like the absence of a spectral parameter, the as- 
sumption to be changed is the order of the underlying scattering problem. 



Let us give more details on the procedure |99|, 100 1 for finding the BT of the 
Boussinesq and Sawada-Kotera equations. 



First example : Boussinesq equation 

Let us consider the Boussinesq (Bq) equation [ 125| , 142 1 

E{u) = uu + e 2 ({u + af + (/3 2 /3)u xx ) = 0, 



(3.50) 



with (a, (3, e) constant. The algorithmic results of the Painleve analysis are : 
p = —2, q = —6, indices -1,4,5,6 compatible 

u T = -2f3 2 X ~ 2 - 2f3 2 S/3 - e-' 2 C 2 /2, X defined by (3.2)-Q (3.51) 

The set of Painleve-Darboux equations reduces to the single equation : 

E 3 = (l/3)/3 2 e 2 S x -C t + CC X = 0, (3.52) 



which is the SM equation for the Bq equation [ 134 1 in the invariant forma-lism. 

(3.53) 



This is a conservation law, which can be parametrised as 

C = ((3e) 2 z x , S = 3z t - (3/2)(Pe) 2 z 2 . 
The compatibility condition of the system ( |3.6| )-( |3T7| ) reads 

3z tt + ((3e) 2 z xxxx + 6(f3e) 2 z t z xx - &(f3e) A z 2 x z xx = 0, 



(3.54) 



which is not the Bq equation but another PDE called modified Bq equation 
[|4[ The elimination of S between ( 3.51 ) and (3.2) yields the Miura trans- 
formation between the Bq and the modified equation 

u T = -(2/3)(3 2 x - 2 - (l/2)s- 2 C 2 + (4/3)p 2 ( X -%. (3.55) 

while the assumption for a DT like ( 3. 29] ) leads to 

u = -(2/3)(3 2 x - 2 - (l/2)e- 2 C 2 - (2/3)(3 2 ( x -%, (3.56) 



which does not coincide with (3.55). We then conclude that a second order 
linear system is not convenient to represent the Lax pair of the Bq equation. 

So, let us assume an underlying scattering problem of the third order for ip 
and the existence of a DT given by the singular part operator 



v T = 2/3 2 LogV> + V, Bq(v T , xx ) = 0, Bq(T4a 



0. 



(3.57) 
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Defining the "second potential Bq" equation 

F(v) = vu + e 2 ((v xx + a) 2 + {(3 2 /3)v xxxx ) = 0, (3.58) 

F(vt) is a second degree polynomial in (Z\, Z 2 ) : 

F(v T ) = C 02 Z 2 + C n Z 1 Z 2 + C W Z 2 + C Q1 Z 2 + C W Z X + C QQ = (3.59) 

which we identify to zero. This provides 

C 02 = 2{{fief - c 2 ) =0, c 2 = (pe) 2 (3.60) 

Cn = -Acd = 0, => d = (3.61) 

C 20 = V xx + a + 2(3 2 a/3 = 0, ^ a = ~3{V XX + a)/(2(3 2 ) (3.62) 

C i = 2{f3e- 1 ac+2{V xx + a) + /3 2 a/3) = 0, => c = /3e (3.63) 

do = 8(/3e) 2 a x /3 + 4e x c = 0, ^e x =/3- 1 eV xxx , (3.64) 

Coo = 2{(3- 2 V xxt + {A/Z)(3sb x + e xx ) = (3.65) 

=> 6 = g(t) - (3/4) (/r 2 ^ + /TV 1 ^). (3.66) 

Finally, the compatibility condition A"o = implies that g(t) is an arbitrary 
constant denoted A. The coefficients a, 6, c, d, e are 

a = -(3/2)/T 2 (J7 + a), c = /?£, d = 0, 

6 = A - (3/4)/r 2 f/ x - (S^)/?-^- 1 ^, e = /T 2 c({/ + a), (3.67) 



i. e. the associated third order Lax pair [141, 142, [56| of the derivative of (3. 
(notation U = V xx ). 



Since d = c x = 0, the BT obtained by eliminating Z 2 between (3.44)-(3.46) 
writes 

Z hxx + 2>Z X Z^ X + Z\ -aZ 1 -b = } (3.68) 
Z ltt + c{Z x Z 1>x + Z\ - aZ l - p- 2 U x - b) = 0, (3.69) 

(Zx, xx ) t - (Z h t) xx = -^- 3 £^(F(V)) X (3.70) 

or equivalently, with U = W x — V xx and Z\ = (w — W)/{2(3 2 ), 

(w - W) xx + 3/T 2 (w - W){{w + W) x + 2a) + /T 4 (w - Wf 
+3(3' 1 £ - 1 (w + W) t -8f3 2 X = 0, (3.71) 
(w + W) xx + [3- 2 {w - W)(w - W) x - p- x s- x {w - W) t = 0, (3.72) 

an extension to A ^ of the bilinear BT of Hirota and Satsuma [|I| [35) . 
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Second example : Sawada-Kotera equation 

In the same way, we can easily find the coefficients of the third order Lax 
pair by processing the fifth order potential equation 

pSK(» =v t + v 5x + 30v x v 3x + 60vl + F(t) = 0, F(t) arbitrary (3.73) 

The algorithmic results of the Painleve analysis are the following : the equa- 
tion (3.73) possesses two families, each with five compatible indices. For the 
"principal" family 

p=—l, q = — 6, uo = 1, indices — 1, 1, 2, 3, 10 compatible (3.74) 
the truncation is 

vt=X~ 1 +vi- (3.75) 



The assumption \ = V'/V'x with ip solution of the second order Lax pair (3.6)- 
(E^tI) generates the Painleve- Backlund equations fj(J 



Eg 



C -AS 2 + 9S XX + Q0Sv i<x - 180vf ia . - 30v l<xxx = (3.76) 
-C x - 2SS X + S xxx + 30S x v ltX - (3.77) 
pSK(m) + {SE A - E Wx )/2 + (5/2)S x (6v hxx - S x ) - 0. (3.78) 



Demanding that v± be another solution of pSK implies v\ tXX = S x /6 and, after 
computation, provides a nongeneric solution. Note that, however, a particular 
solution of the truncation is 132 



S xx +4S 2 -C = 0, v hx = 5/3, KK(»x) = 0, 

which defines a Miura transformation between SK and KK equations. 
As in the preceding example, the hypothesis of the DT 



(3.79) 



u = (LogV0x+V, (3.80) 

with V another solution of pSK and ip a solution of the third order linear 
system (3.37)-(L3S), makes pSK(v)— pSK(F) a second degree polynomial in 
(Z±, Z-i) like (3.59). The six determining equations C; m = 0, added to the three 
compatibility conditions ( 3.40| )-( |3.43D , have the unique solution depending on 
an arbitrary constant A 

a = -6V X , b = A, c = 9A - 18V XX , d = -36V 2 + 6V 3x , e x = 36\V XX , (3.81) 



a result which coincides wi th the Lax pair ( |2.144| ). 

T he a:— part of the BT ( 2.148 ) is obtained by eliminati ng Z 2 between ( 3.44 ) 
and ( 3.45 ), then substituting Z\ = v — V as results from ( 3.80 ). 
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Some results for Kaup-Kupershmidt equation 



In the case of potential KK equation, the hypothesis of the di fferen tial op- 
erator T> — d x for the DT, associated to the linear system ( 3.37)— (3.38), yields 
neither the Lax pair (2.145) nor the BT ( 2.149] ). This problem has been recently 
solved [102 1 by remarking that in his classification of second order first degree 
nonlinear ODEs possessing the Painleve property, Gambier |54j] mentions that 
the following equations : 



(G.5) : 
(G.25) 



A' 



Y 2 Y 2 ^ XX — 7^2, X 



Yf + WY X - A = 
1 



2 * -~r» ■ 4 

are linearisable into third order equations 



•JffYi* + 7*7 + §UY4 - 2\Y 2 = 0, 



(G.5) : Yi-. 
(G.25) : Y 2 

ipxxx + 6Uip x + (3U X - X)ip = 



= ip x /i/>, i/j xxx + 6Uip x - Xip = 
- 1 = X- 1 [(tP x /^) x + (l/2)^ x /i;) 2 



■3J7], 



(3.82) 
(3.83) 

(3.84) 
(3.85) 



corresponding to the scattering problem of, respectively, the SK equation for U 
and the KK equation for U. It can then be shown that the DTs 



Yi = w - W, with SK(w x ) = SK(W X ) = 
Y 2 = 2(w - W), with KK(w x ) = KK(W X ) = 



(3.86) 
(3.87) 



leading to the BTs (2.148) and (2.149), can be found by singularity analysis. 



3.2.4 Two-singular manifold method 

For equations with two families of movable poles with opposite residues, the 
truncation procedure which considers only one family of singularities does not 
yield the auto-BT. An extension of the Weiss method consists of considering 
two distinct functions -01 , "02 , assuming now a DT of the form 



u T - U = X>Log0i - VLogfa, E(u T ) = E{U) = (3.88) 

and, assuming that Y — ip\ /ip2 satisfies the most general Riccati system 

Y x = Ra + R^ + ^Y 2 (3.89) 

Y t = S a + S 1 Y + S 2 Y 2 (3.90) 

Y xt -Y tx = X Q + X 1 Y + X 2 Y 2 , (3.91) 

Xq = Ro : t — Sq :X + RiSq — RqSi, (3.92) 

X 1 ee R ht - S hx + 2(R 2 S - R S 2 ), (3.93) 

X 2 ee R 2t — S 2 . x + R 2 S\ — RiS 2 , (3.94) 
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eliminating derivatives of Y and identifying Et to the null polynomial in Y. 

The determining equations so generated for the six unknowns (Ri, Si) must 
have a solution such that each Ri is a linear function of U and an arbitrary 
constant A and such that at least one of the three cross-derivative conditions 
Xi = is identical to the original equation for U. In such a case, one has 
found the DT and the Lax pair, i.e. the BT. In the case of two simple poles with 
constant opposite residues ±ti an d opposite singular part operators ±u d x , the 
truncation 

u T = u (hogY) x + U (3.95) 



becomes, by elimination of Y x from ( 3.89 ) 



ut = uo{R Y- 1 +R 1 + R 2 Y) + U. (3.96) 

This represents an extension of the Weiss truncation to the positive powers of 
Y. A similar extension was previously made ]109[ in the variable \ = ip/rfx 
for obtaining particular solutions of nonlinear PDEs. The "two-singular mani- 
fold" method is successful for finding the auto-BT of MKdV and sine-Gordon 



equations 1 101 ] but only partially for the NLS equation. Before detailing this 



result, let us first reproduce in Table 3.1 the algorithmic results of the Painleve 



analysis for the four equations belonging to the AKNS scheme; these include the 
SM equation associated with these well known NLPDEs which pass the Painleve 
test. 



Table 3.1: Algor ithmic results of the Painleve anal ysis. T he integers (p, q) are 
defined in ( 3.26 ); for sG, the polynomial PDE is ( |3.108 ). Next column lists 
the indices, excepted —1. Column "PD equations" lists the subscripts of the 
non identically zero Painleve-Darboux equations; in the sG and NLS cases, they 
depend on the arbitrary coefficients introduced at the index 2 (sG) and (NLS). 



Name 


P 


q 


Indices 


PD eq. 


Singular 
manifold 
equation 


KdV 


-2 


-5 


4,6 


3,5 


S - C = 6A 


MKdV 


-1 


-4 


3,4 


2 


s-c = o 


sG 


-2 


-6 


2 


3,4,5,6 


S + C i C xx 
-C- 2 Cl/2 
+2A = 


NLS 


(-1,-1) 


(-3,-3) 


0,3,4 


2,2,3 


Ct + iCC x S x 
+8\C X = 



Let us now use the information contained in the SME. 
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Modified Korteweg-de Vries 



The equation ([2.93|) has two families u ~ ±a% 1 , denoted u ~ ax 1 since a 



is defined by its square. The truncated expansion of a family is 

u T = ax" 1 - 
The SME S — C = is parametrised as 

S = 2v, C = 2v, KdV(v) = 0, 



(3.97) 



(3.98) 



and the precise relation between u and v (Miura transformation) is obtained by 
eliminating x between ( 3.97 ) and ( |3.2| ) 

(u T /a) x + (u T /a) 2 = -v. (3.99) 

In fact there are two such Miura transformations, one for each sign of a, i.e. one 
for each family. 

Let us first obtain the Darboux transformation for MKdV from that of KdV 
and show that it involves two SMs. The Darboux transformation for KdV has 
been obtained in section 3.2.2, eq. ( 3.36| ). The two Miura transformations (3.99) 
and the parametrisation (3.98|) imply : 



Si _ ^u T y | (u T ^ 



S2 

2 



V a 



(-) 



2(Log-0i) a 
2(Log^2) a 



(3.100) 
(3.101) 



and the elimination of the nonlinear terms leads to 

ut.x = a(hog(tpi/ip 2 ))xx + U x , 



(3.102) 

which after one integration yields the Darboux transformation for MKdV 

u T = a(Log{i/fi/ip2)) x + U. (3.103) 
With this DT, the Lax pair is obtained as explained in the introduction 



of this section. Setting Y — tpi/4>2 and taking account of ( p.89[ )-( p.90| ), ev- 



ery derivative of Y can be replaced by a polynomial in Y . Consequently, the 
Darboux transformation (3.103) becomes identical to 



u T = aiRoY- 1 +Ri + R 2 Y) 
and one must identify to zero the polynomial in Y 



U 



(3.104) 



E(u T ) = E t = Y^ E J Yj ~ 

3=0 



(3.105) 
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Among the nine Painlevc-Darboux equations Ej = 0, only four (J — 1,2,6,7) 
are not identically zero. Their resolution, as detailed in [ 101 1 , yields the following 
parametric representation of the six unknowns Ri , Si in which Ri is linear in U 
and the spectral parameter A 

R a = A, R 2 = -A, Ri = -2U/a (3.106) 

So = -AX 2 + 2{U/a) 2 -2U x /a, Si = 8X 2 U/a - A(U/a) 3 + 2U xx /a 

S 2 = -4A 2 + 2{U/af + 2U x /a. (3.107) 

This soluti on as so ciate d with the Riccati equations ( 3.89| )-( 3.90 ) reproduces the 
equations ( 2.98|) -( ^99 ) for the pseudopotential of MKdV . 

Sine-Gordon 

The sine- Gordon equation ( p. 118 ), invariant by parity on u, is first trans- 
formed into a polynomial equation for v = e m , invariant under v — > l/v. 



PsG(u) = 2vv xt - 2v x v t - v 3 + v = 0, v — e 1 



This PDE has two families of movable singularities v 



(3.108) 
-4Cix~ 2 and 



-4C2X ■ The truncation equations have the following general 



solution [133, E9J]. For the first family ((v, S, C, ip) are subscripted with 1) 



Si = -v hxx /vi+vlJ(2v 2 1 )-2X-- 

Ci = -«i/(4A) - -e m /X 

vi = -AfLogfaU + Vx, PsG(Vi) 



-iu xx + u 2 /2 - 2X 



For the second family e 



v 2 



S-2 = -v 2 , xx /v 2 +v 2 1 x /(2v 2 1 )-2X = iu xx + ul/2-2X 

C 2 - -«a/(4A) = -e- iu /X 

v 2 = -4(Log^ 2 ) xt + V 2 , PsG(V 2 ) = 0. 



(3.109) 
(3.110) 
(3.111) 



(3.112) 
(3.113) 
(3.114) 



If one considers only one of these two equivalent SMs, the Schwarzian Si,i — 
1 or 2, does depend on an arbitrary constant A but it has two drawbacks : it 
is not invariant under parity on u, it is not linear in the physical field u as 
requested for the Lax pair (3.6)-([3j]) to be a "good" one. 

Since v± — v 2 — 2isinu, the difference of ( [3.111 ) and ( [3.114 ) reads 



sinu = 2i(Log(^i/V> 2 ))a:t +sinJ7, sG(U) = 0, 
i.e., from the definition of the equation 

u xt = 2i(Log(ipi/ip 2 )) xt + U xt . 



(3.115) 



(3.116) 
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Integrating twice, we finally obtain the Darboux transformation of sG 

u = 2iLog(^i /ip 2 ) + U, (3.117) 

defined in terms of both families. For the solution of the polynomial PDE 
PsG(v) = associated to the sG equation by v = e m , the DT is : 

v = VY-' 2 , Y = ViM, PsG(V) = 0, (3.118) 

and one must identify 

E(v) — V 2 y ' E t Y % ~^ (3.119) 



4 
i=0 



to the null polynomial in Y. Among the five Painleve-Darboux equations, the 
equation E% is functionally dependent on (Eg,Ei), a consequence of the com- 
patibility of the index 2. Their resolution yields the Riccati pseudopotcntial 
( |2.122D -( |2.123|) 



Y x = A(l - Y 2 ) + iU x Y, (3.120) 

Y t = ((l-r 2 )cosC/ + i(l + r 2 )sinC/)/(4A), (3.121) 

(Y xt -Y tx )/Y = sG(U), (3.122) 

where the x-part is now linear in the spectral parameter A and the field U 
associated with the DT fl3.117| ). 



3.2.5 Weiss method plus homography 



Pickering [11C] remarks three drawbacks in the previous method 



(i) any explicit relationship is given between the variable Y = ipi/ip2 and the 

variable x of the invariant Painleve analysis while the one between \ & n d 
if is well defined by the homographic transformation pgj 

X= ~ 775 y (3-123) 

(ii) the result of the previous truncation for MKdV does not reveal any rela- 
tionship between the MKdV and KdV equations as it would be, 

(iii) the knowledge of the DT is required in advance. 

He notices that for finding the BT of MKdV and sine- Gordon equations it is 
sufficient to consider a Riccati system constructed from the nonlinearisation of 
the following second order scalar linear system 

r) xx = 2Ar) x + Br) (3.124) 
Vt = -Cr h +( f D da/) V (3-125) 
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by the transformation Z 1 = i] x /r]. The corresponding nonlinear system 



Z X = 1-2A- BZ 2 

Z t = -C+ (C x + 2AC)Z - (D— BC)Z 2 



(3.126) 
(3.127) 



depends on four functions A, B, C, D in place on six like the system (|3~89|)-(|3~90|) 
considered in the 2-SM method. Its compatibility condition is 



Z x t — Z tx — X\Z + X2Z 2 — 0, 

X 1 = 2 (D - {A t + (AC) X + C xx /2)) = 

X 2 = D x - B t - 2BC X - B X C - 2AD = 



(3.128) 
(3.129) 
(3.130) 

The solution 77 of the linear ODE (|3.124| ) is related to the solution ip of ( |3.6| ) by 
the gauge transformation 



(3.131) 



with 



S = -2{B + A 2 - A x ) (3.132) 
Then, computing the a;— derivative of Log 77, one gets the transformation 

Z~ 1 =x~ 1 +A (3.133) 

which means that the new expansion variable Z is related to x by an homo- 
graphic transformation (as suggested in |9{J, formula (17)) su ch tha t in the 
neighbourhood of % = 0, one has Z ~ \- Then the system ( 3.126 )— ( 3.127 ) 
combined with the truncation in Z 



ut = ad x Log Z + [/, 



(3.134) 



(U function of (x,t) and a constant) such that u T ~ ax 1 as X 0, extends 
the Weiss truncation to positive powers of Z. Solving the Painleve-Darboux 
equations associated with E(ut) = Pickering obtains for the MKdV equation 
( [2.9S ) the following results : 



A = U/a, B = A 2 , C = 2(U x /a - (U/a) 2 + 2A 2 ) 

D = 4X 2 U x /a, A = A(i) arbitrary integration function 

and the compatibility conditions ( |3. 129 )— ( 3. 13C ) yield 

Xi = -(2/a)(U t + U xxx -2a- 2 (U 3 ) x ) 
X 2 = -(A 2 ) t -0 

From (^.132[) and (|3.135|) one gets 

S , = 2(C/ :r /a-(C//a) 2 -A 2 ) 







S - C + 6A 2 



0. 



(3.135) 
(3.136) 



(3.137) 
(3.138) 



(3.139) 
(3.140) 
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the latter equation being the SM equation of the KdV equation. 

Let us remark that the expressions of S and C in function of the solution U 
of the MKdV equation and the constant parameter A coincide with the relation 
( p.lll ) obtained previously. The expression (2.10S) for Y with a = U/a implies 
the identification Z = \~ 1 Y. 

For the sine-Gordon equation ( 2.11S| ), considering the truncation 



Ut = 2i Log Z + U 
such that ut ~ 2i Logx as x ~~ * 0, Pickering obtains 



A = — (i/2)U x , B = X , D = —[i/2)smU, A arbitrary constant, 
C = -X~ 2 e iU /4, S = -iU xx + (l/2)U 2 - 2A 2 , sG(f/) = 0, 
5 + C xx /C - {1/2){C X /Cf + 2A 2 = 



(3.141) 



(3.142) 
(3.143) 
(3.144) 



Again the expressions for S, C in function of U and A and the equation (3.144) 
coincide with the relations ( 2.12S| ) and (2.13C). The identification Z = A _1 Y" is 
also easy to find taking account that, for sG, a — —iu x /2. 



3.2.6 Weiss method plus involutions 

The AKNS system |l39j § 



lU t + p r U xx 



+ q r u 2 v = 0, = —iVf + p r v xx + q r uv 2 = (3.145) 



has the BT ||, H |7| |§ (a 2 = -2p r /q r , R 2 = {u + U) (v + V) /a 2 - (A — /.i) 2 ) 



(« + U) x 
(v + V) x 
-ip~ x (u + U)t 

-ip-\v + V) t 
M 



= -(u - U)R - i(X + n){u + U) 
= -(v-V)R + i{\ + n)(v + V) 
= (it - U) X R +{u + U)M + i(X + fi)(u + U) x 
= (v - V) X R +(v + V)M - i(X + fj,){v + V) x 
= {uv + UV)/a 2 



(3.146) 



with A, fi arbitrary complex constants. Galilean invariance (x, t, u, v) - 

■M = oJra|. 



2p r ct, t, e t( - cx -P- c Qu, e -i(^-Prc t) w ) a n ows to choose c = A 



The above BT can not be found neither by the one-SM method [134], n or b y 
the two-SM method [101|, nor by the one-SM method plus homography |11C|. 
The challenge of the Painleve approach to find this BT by singularity analysis 
only is solved in 1 35 as follows. 



As the one-SM method only provides some partial result T(x, u, X) for the 
truncation, one then considers all transformations on u conserving the equation 



E(u) = in order to unco ver a s econd solution U, see Table 3.2 
For the AKNS system (3.145), the one-family truncation 



u Q x + Ml, v = VoX + Vl 



(3.147) 



39 



Table 3.2: Transformations of the dependent variable(s) conserving the equa- 
tion^), for the AKNS group PDEs (complex conjugation, phase shift, parity). 



PDE 


Transformation(s) 


AKNS system 


(u, v, i) — > (v, u, —i); Vfc : (u, v) — > (ku,v/k) 


Sine-Gordon 


u — ► — u 


MKdV 


u — > — u 


KdV 


none 



which has the general solution |134, |101[| (A arbitrary complex constant 



u = aix- 1 - f x /(2f) - iX)f 
v = a( X - 1 +f x /(2f)+i\)/f 
f x /f = -2iX-(u/a)f- 1 + (v/a)f 

*JV Vt// = 2uv/a 2 + 4A 2 + {u x - 2iXu)/(af) + (v x + 2i\v)f/a 
(fxt-ft*)/f=tf- 1 E<»+fEW)/a 



(3.148) 
(3.149) 
(3.150) 
(3.151) 
(3.152) 



fails to introduce a second solution {U, V), sec details in appendix C of Rcf. [101 1 . 
This is done by applying the two point transformations of Table 3.2 to the above 
truncation T x : 



(3.153) 



Tx : 


Xi 


u 


V 


* f 


X 


(identity) 


T 2 : 


X2 


V 


u 


-i 9 


fJ> 


(conjugation) 


T 3 : 


X3 


kU 


k- x v 


* f 


X' 


(phase shift) 


T 4 : 


Xa 


k~ Y V 


kU 


-i 9 


K' 


(both) 



These transformations act on (u,v,f,X) like in Chen p3| . This is equivalent 
to successively process the four families of the AKNS system by the one-SM 
method. In order that (u, v) and (kU, V/k) be distinct, one must have A' = 

The four sets ( |3.148| )-( |3.149| ) define a system of eight equations in the eight 
unknowns (Xi^'X^ 1X3^1X4 ,u,v,kU,V/k). This system is linear with deter- 
minant fg — l/(fg) and it provides the DT straightforwardly (with the nonre- 
strictive choice k = — 1) : 



u-U = 2a[d x Log(g-l/f)-i(X 

v-V = 2a[d x Log(f-l/g)+i(X- 

u + U = 2m(A -//)/(<? -1//) 

v + V = 2ia(X-fi)/(f-l/g) 



■A»)]/(/- 



V/) 
1/9) 



(3.154) 
(3.155) 
(3.156) 
(3.157) 
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(to stick to our definition, the DT is made of two equations, either ( |3.154 ) 
( |3. 155 ) or (3.156)— (3.157)). The nonconstant factor of the logarithmic deriva- 
tives is similar to that of (P3), (P5), (P6), see section 7.1 in the Conte lecture 
while A + /i is choosen as a real constant. 

The Lax pair in its Ric cati fo r m is ma de of the four equations resulting from 
the action of T 3 and T4 on ( 3.150 )— ( 3.151 ). The BT is made of the four equations 
resulting from the elimination of the pscudopotcntials (/, g) between the six 
equations defining the DT and the Lax pair, and th es e are precisely ( 3.146 ). 
This elimination is quite easy since equations ( 3.156| )— (3.157) are algebraic in 
(/,<?) : 

/ = ia{\ -n + R)/(v + V), g = ia(X -fi + R)/(u + U). (3.158) 

Remark. The system ]23| of t wo equations for (/, g) , obtained by elim- 
inating (u,v, U, V) between ( |3.154 )— ( 3.157 ) and the PDE, is invariant under 
(A,/i) — * (fx, A). The elimination of g between this system provides the Brocr- 
Kaup equation for w = —i Log /, a result also obtainable by the Weiss trunca- 
tion Jl01| 



p r 1 w tt + <iw x w xt + 2w t w xx + p r (6wlw xx + w xxxx ) = 0. (3.159) 



3.2.7 Reductions of the DT of AKNS system 

The x— part of the AKNS spectral problem admits the three reductions v = 
u, v = ±m, v — 1, and the DT, obtained only from the x— part, must admit them. 
This is indeed the case : equations (3.154 )— ( 3. 157) admit the two reductions 
(v, V,g, /i) = (u,U,f,\),(eu,eU 7 sf,—\),£ 2 = 1, and one must add the case 
g = e/f when the determinant vanishes. Table 3.3 summarises these reductions 
and the homographic link between / and the x of the invariant analysis. 
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Table 3.3: Reductions of the Darboux transformation of the AKNS system. 



PDE 


v,g,n 


X" 1 


(u-U)/a 


(u + U)/a 


NLS 


u,f,X 




(|3.154|) 


4(Im A)/(l// - /) 


sG 

MKdV 


eu,ef, -A 
e 2 = e 


A eU 

Y ~ " 4o ' 
Y= ef-1 
e.f + 1 


(4/e)£ = 

£f 2 -l 


e y 

= 4*A— ^— - 
e/ 2 - 1 


KdV 


M//,-A 


f-iX 


-2(X" 1 ), 


2{f - 2iXf) 
= 2( x - 2 + A 2 ) 



3.3 Cosgrove classification for semilinear PDEs 
of second order 

In two papers |37], [38| , Cosgrove classifies two cases of Painleve type semilin- 
ear PDEs of second order. The necessary conditions for the Painleve property 
which he establishes combine the criteria of Painleve and Gambier for ODEs 
and the WTC ones for PDEs. 

For hyperbolic PDEs in two independent variables of the type 

u xt = F(x,t, u,u x ,u t ) (3.160) 
the classes of equivalence are defined by the H-transformation 

a(x,t)u + f3(x,t) . i el \ 

U = — -, : —, (3.l0l 

j(x,t)u + S(x,t) y ' 

whereu = u(x, T), x = X(x), t = T(t), aS - (3j ^ 0. (3.162) 

and the necessary conditions are 

1. the dependence in u x ,ut of F must be of the form 

u x t = A(x, t, u)u x ut + B(x, t, u)u x + C(x, t, u)ut + D(x, t, u), (3.163) 

2. as a function of u, the term A is the sum of at most three simple poles, 
at locations set to u = 0, l,H(x,t) (H arbitrary function of x, t) while 
B,C,D cannot grow faster than, respectively, u,u,u 3 , 

3. the equation must pass the WTCK Painleve test, i.e. all Fuchs indices 
are distinct integer and all positive indices are compatible in order to 
guarantee the existence of local Laurent series in the Kruskal variable 
ip(x, t)=x± fx(t) or ip(x, t) =t± f 2 (x). 
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At the end, he obtains 22 canonical equations, reducible to 

u xt = sin it (sine-Gordon) or u xt = ae 2u + be~ u (Tzitzeica) (3.164) 

or linearisable by the singular part transformation. 

For parabolic PDEs the class of equivalence is a little bit larger than 
in the previous case in the sense that the new independent variables in the 
transformation ( [3.161 ) may be related to (x,t) as 

x = X(x,t), t = T(t) (3.165) 

The successive necessary conditions for having the Painleve property yield the 
following results 

1. in two independent variables, the sole equation is 

u t + u xx + 2uu x — F(x,t), F(x,t) arbitrary function (3.166) 
i.e. the Forsyth-Burgers equation linearisable into the heat equation, 

2. in more than two independent variables, only trivial soliton equation are 
obtained, i.e. nonlinear PDEs related to linear ones. 



3.4 PDEs with variable coefficients 

The Painleve test can be applied to nonlinear PDEs with variable coefficients 
in order to determine the conditions under which the equation might be intc- 
grable. The sufficient part of the analysis entails the determination of the DT 
and the Lax pair, as well as the transformation which could relate the equation 
to its autonomous integrable counterpart. Many authors have considered the 
generalised KdV and NLS equations with variable coefficients due to their inter- 
est in many physical systems. The results of Brugarino |p^| for the generalised 
variable coefficient KdV equation (VCKdV) 

u t + a{t)u + (b(x, t)u) x + c(t)uu x + d(t)u xxx + e(x, t) = 0, (3.167) 

are 

(i) the equation passes the Painleve test under the condition 

b t + (a — Lc)b + bb x + db xxx = 2ah + hL(d/c 2 ) + bl + ce 

+x (2a 2 + aL{d 3 /c 4 ) + a' + L(d/c)L(d/c 2 ) + (L(d/c))') (3.168) 

with L = (d/dt) Log and h(t) arbitrary, 

(ii) the solution of the Weiss truncation yields the DT and the Lax pair, 
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(iii) with the transformation from (u,x,t) to (0,£,i) 
u = ((a + L{d/c))i + g-b + 0)/c, £ = x- 



g{T)AT (3.169) 



and the condition ( 3.168 ) , he gets the equation 

Qt + (2a + L(d/c 2 )) 6 + £ (a + L(d/c)) 6 e + 69^ + d6 5ee = 0, (3.170) 

equivalent to |^l|, [l7], |66| the KdV equat ion with constant coefficients. In 
case a = b = e = 0, the condition ( [3.168] ) becomes simply 



d = c (kx J c(T)dT + K^j , K 1 ,K 2 arbitrary constants, (3.171 



i.e. the one given by Joshi p7| when performing the Painleve test on 
this particular VCKdV equation. The same relation was obtained by 
Winternitz and Gazeau [ 137 1 using the symmetry group, 



(iiii) several equations of physical interest which satisfy ( 3.168 ) are presented. 

Gagnon and Winternitz p3| have analysed from the point of view of sym- 
metries a variable coefficient nonlinear Schrodinger (VCNLS) equation 



iu t + f(x, t)u xx + g{x, t)u\u\ + h(x, t)u = 



(3.172) 



involving the three complex functions /, g and h of the variables x, t. The sym- 
metry group is shown to be five-dimensional iff the equation ( 3.172| ) is equivalent 
to NLS itself or to CGL3, and at most four-dimensional in all other cases. In this 
framework, they give the allowed transformations of (3.172), in case / = 1 + 4/2, 
to the equation with constant coefficients 

iu~ t + (1 + if2)Uxx + (h + «.92)"|u| 2 + {hi + ih 2 )u = (3.173) 

which are, in the special case f% = (/, J, K, L arbitrary functions of t), 



9 = (§1 + ifa)TI 2 

Reh=(K + AK 2 )x 2 + (L + AKL)x + j + L 2 + hit 
Im/i = -ir 1 - 2K + h 2 f 

t = T, f = T e~ 8 f Kdt 1 x = 
u = u(x,t)Ie t( - Kx2+Lx+J K 



Tx + £„ £ = — 2v TL 



(3.174) 
(3.175) 
(3.176) 

(3.177) 
(3.178) 



One example given by the authors leading to the NLS equation {g 2 = h\ — 
h 2 = 0) corresponds to the following choice of the arbitrary functions involved 
in the transformation 



IT 1 = -2K, J = L = 0, K + AK 2 = Kq/4, K constant 



(3.179) 
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(3.180) 



and yields the VCNLS equation 

iut + u xx + cjiTo sech(\/ Kot)u\u\ 2 + {Kq/ 4)x 2 u = 0. 

This equation is related to NLS by the change of variables 

x = Tj%xsech(T,/K^t) + x , t = T tanh( v / ^o*)/v / ^o (3.181) 
u{x, t) = u(x, t) sech 1/2 ( y/Ktf) e^/^v^ 2 tanh(v^t) ( 3182 ) 

Considering the other choice : TI~ 2 = 1, one obtains 



ilk + u xx + giu\u\ + (AK 2 + K)x z + (L + AKL)x + L + J + 2iK )u = 



equivalent to NLS by the transformation 



i = To / At' e 8 $ Kds . x = x\/T e 

£ = Co - 2^ / &t'L(t')e- A r K(s)As 
u{x,t) = u(i 1 t)e- i r Kdt 'e l(Kx2+Lx+J '> 



(3.183) 



4 J* Kdt' 



(3.184) 



For 2f(t) = -/?(*) /2, Equ. ( [3.183; ) and the transformation ( |3.184| ) coincide with 
the results of Clarkson [E8J in the analysis of the PDE 



iu t + u xx — 2u\u\ = a(x, t)u + b(x, t) 



(3.185) 



which passes the Painleve test iff there exist functions (f3(t),ai(t),ao(t)) such 
that 

a = i/3 + x 2 [/3'/2 - (3 2 } + xa x + a , b = 0. 



As noticed by Clarkson, this result proves that the equation 

iut + u xx — 2u\u\ 2 = f3x 2 u, [3 constant, (3 ^ 0, 
which does not satisfy the condition ( 3.186| ), is not integrable while 

iu t + u xx — 2u|w| 2 = (in — k 2 x 2 )u, k real, 
can be transformed into NLS and hence is integrable. 



(3.186) 
(3.187) 

(3.188) 
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Chapter 4 

Partially and non integrable 
equations 

Let us consider a polynomial PDE (depending on a parameter /i) which does 
not pass the Painleve test but possesses a singularity structure compatible with 
singlevalued solutions. The reason for the failure of the Painleve test are the 
following ones : 

(i) the Fuchs indices are all integers (possibly for particular values of fi) but 

some of them do not satisfy the compatibility condition, therefore the 
existence of a Laurent series is submitted to conditions on the SM, 

(ii) whatever be jj,, some indices are irrational and thus unable to generate com- 

patibility conditions, therefore the number of arbitrary functions which 
can be introduced in the Laurent series is lower than the order of the 
equation; this generally characterises equations with chaotic behaviour. 

In the first case, we classify the equation as being partially integrable while in the 
second case as being nonintegrable. We provide methods for finding particular 
closed form solutions and illustrate them on some examples. In each case one 
looks for solutions related by a rational transformation to the general solutions of 
first order nonlinear ODEs like Riccati, Weierstrass or Jacobi possessing the PP, 
or solutions related by a nonlinear transformation like the logarithmic derivative 
to a second or third order linear system with constant coefficients. 

4.1 Partially integrable equations 
4.1.1 KPP equation 
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The Kolmogorov, Petrovskii and Piskunov equation (KPP) |7l], 103 1 



E — u t — u xx + ^p( u ~ e i)(u — S2)(u — 63) = 0, ej distinct, (4-1) 
possesses two opposite families of singularities 

u = dx' 1 + si/3 - dC/6 + O(x), indices : - 1, 4, V = dd x (4.2) 
s\ = ei + e 2 + e3, s 2 = e 2 e3 + + eie 2 , S3 = eie 2 e3. (4.3) 

with the condition on the SM coming from the index 4 

3 

Q 4 = C[-3(C t +CC s ) + *[[{C + 8i- 3e fc )] = (4.4) 

k=l 

Denoting (j, Z, m) any permutation of (1, 2, 3) and fcj = (3e^ — Si)/(3d), the one- 
family truncation u = dx^ 1 + si/3 — dC/6 yields the moving front solution : 

u = - — — + rf— tanh — (x — ct — in) (4.5) 
2 2 2 

fc 2 = {ki-k m f, c=-3(fc ; + fc m )/2. (4.6) 

The two-family truncation like for MKdV or sG yields the stationary pulse 

. e 2 - e 3 e 2 - e 3 

u = ei H = — sechi -^(x — xo), 2ei — e 2 — ez — 0. (4.7) 

\/2 dV2 

A third very interesting solution representing the collision of two fronts j70) with 
different velocities can be easily found |33| with the assumption 

u = y + u d x Logip, (4.8) 

where ip is the general solution of a linear system with constant coefficients 

3 

ip = C ™ exp [k„(x + b 2 t) + klht] , C„ arbitrary dC 2 C 3 ^ (4.9) 

n=l 

ipxxx - a\i>x - a.2i> = 0, ip t - biip xx - b 2 ip x = (4-10) 

with the following values of the coefficients and the constant k n 

ai = (s 2 - 3s 2 )/(3d 2 ), a 2 = [2s\ - 9s lS2 + 27s 3 )/(27d 3 ), 

h = -3, b 2 = 0, k n = (3e„ - Sl )/(3d). (4.11) 

4.2 Nonintegrable equations 
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4.2.1 Kuramoto-Sivashinsky equation 

The Kuramoto-Sivashinsky (KS) equation 

E = m + uu x + u xx + u xxxx = 0, (4-12) 

describes, for instance, the fluctuation of the position of a flame front, or the 
motion of a fluid going down a vertical wall, or a spatially uniform oscillating 
chemical reaction in a homogeneous medium. For a review, see ref . ]89[ . 
This equation possesses only one family of singularities |5lJ 

u = 120x~ 3 + 60(5 + l/19)x _1 + (C - 15S X ) + 0( X ), (4.13) 
indices: - 1, 6, (13 ± iVn)/2, V = 60d x 3 + (60/19)5* (4.14) 

Due to the existence of the two complex irrational Fuchs indices the Laurent 
series depends only on two arbitrary functions (uq and the arbitrary function in 
the expansion variable \) whatever be the SM. The WTC truncation 

u = 120x~ 3 + 60(5 + l/19)x -1 + (C - 15S X ) (4.15) 

generates three determining equations E4 = 0, E$ = 0, £7 = whose general 
solution is C = arbitrary c, 5 = —11/38,1/38. This corresponds to the well- 
known travelling wave solution of Kuramoto and Tsuzuki [f79[ only existing for 
two values k 2 = -1/19 or 11/19 : 

30 \ k k 

— k - 15/fc 3 tanh -(£ - x ) + 15fc 3 tanh 3 -(£ - x ) (4.16) 
19 J 2 2 

where £ = x — ct and c, xo are arbitrary constants. 

This solution can also be retrieved with the assumption that u = c + T> Log ip 
with ip the general solution of a linear system with constant coefficients 

- (fc 2 /4)V> = , ik - ap x = 0. (4.17) 



Let us remark that the reduction u{x,t) — > c+ U(x — ct) of the PDE ( 4.12 ) 
yields the nonintegrable ODE 

U'" + U' + U 2 /2 + K = 0, K arbitrary (4.18) 

for which we have found, in the c ase K = — 450/c 2 /19 2 , a Riccati sub-equation 
linearisable into the system (4.17). The challenge not yet solved is to find for 
every K a closed form particular solution depending on one arbitrary constant. 

Let us also mention the interesting work of Porubov [111, 112 who has found 
for a large class of nonlinear PDEs like 

T)t + aiVx + a 2 r/Vx + a 3 (r/Vx)x + a^xx + a 5 r] xxx + a 6 r)xxxx = (4.19) 
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and for some particular values of the constant parameters {<2i}, travelling wave 
solutions in terms of Weierstrass elliptic functions or its logarithmic derivative. 



For the KS equation with an additional dispersive term 

E = u t + uu x + u xx + bu xxx + u xxxx = 0, (4.20) 



which belongs to the previous class ( |4.19| ), Kudryashov |75j has found a partic- 



ular solution depending on three arbitrary constants in terms of the Weierstrass 
elliptic function and its derivative in the single case b 2 = 16. An easy way to 
find it is to assume the existence of a solution of the form 

u = a Q + a 2 p{x - ct- £0,32,53) + a 3 p'(x - ct - 2:0,32, 93), (4-21) 

compatible with the singularity structure of its Laurent series expansion. En- 
forcing in the LHS E the conditions p' 2 = 4p 3 — 32 p — 33 and p" = 6p 2 — (32/2), 
one identifies to zero a polynomial of two variables (p, p') of degree one in p'. 
This similarly generates six equations in the five unknowns a,j , 32 , 33 and yields 
a nondegenerate elliptic solution only if b 2 = 16 and 

ao = c — 4/6, ci2 = —156, 03 = —60, 32 = 1/12, (33,0,2:0) arbitrary. (4.22) 



4.2.2 Complex Ginzburg-Landau equation CGL3 

The one-dimensional cubic complex Ginzburg-Landau (CGL3) equation 

E = iu t + pu xx + q\u\ 2 u — i"fu = 0, pq 7^ 0, (u,p, q) £ C, 7 € 1Z, (4.23) 

with p, q, 7 constants, describes pattern formation and coherent structures in 
many different domains : Taylor-Couette flows between coaxial rotating cylin- 
ders, wave propagation in optical fibers and chemical reactions. For a review 
sec Q. This PDE is physically strongly connected to the KS equation and it 
also possesses two irrational complex conjugate indices. As is the case for the 
AKNS system, CGL3 equation possesses four families of singularities 

u = A QX - 1+ia (l + A lX + 0{ X 2 ) ) u^Box' ^jl + B 1X + 0( x 2 )), (4.24) 

A B = 3\p 2 \a/D u a = D r ± ^ D 2 r + 8D 2 /9, A = p r<?i - Pi q r , (4.25) 

A = (PrQr +PiQi), Fuchs indices : - 1,0, 7/2 ± \J 1 - 24a 2 /2. (4.26) 

The important information we get from the singularity analysis is that neither 
(u,u), nor (|it|, argil) nor (Re u,Imtt) have a simple singularity structure. In 
this framework, better variables are (Z, 6) defined as 

u = Ze ie , 9 = a Logip + e Q , </>M = X, (4.27) 
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where 9o is an arbitrary function representing the index 0. Then Z and grad# 
behave like simple poles, so that the usual methods are applicable. With tp the 
general solution of a second order or third order system and with a one-family 
or two-family truncation for [Z, 9) one finds |33| all the known closed form 
solutions of this equation. Among th em, a very interesting solution representing 
a "collision of two shocks" [104, 105 



A k 



sinh kx/2 



2 coshfcx/2 + e- 3 7(t-*o)/2 

k 2 = -2 7 / Pi , Pr = 0, 



i[a Log(l+e 3 "' (t " t o' /2 cosh kx/2)] 



(4.28) 



is easily found by assuming that ip satisfies a third order linear system. For 
q r = (hence a = 0) this solution degenerates to the "collision of two fronts" 
solution of KPP previously considered. The Table 4.1 summarises the known 



solutions of CGL3 and their degeneracies to NLS and KPP equations. An open 
problem is to find the solution of CGL3 which degenerates for Pi = qi = J — 
to the bright soliton ( |2 27| ) of NLS. 



Table 4.1: Degeneracies of the known solutions of CGL3. 



CGL3 


NLS ( Pi = qi = 7 = 0) 


KPP (p r = q r = 0) 


propagating hole |9, 33] 


dark soliton (2.28) 




shock or front 1 105 , 33J 




front (4.5) 


pulse or solitary wave ]fl08|, |33| 




stationary pulse (|4.7|) 


collision of two shocks (p r = 0) 
(4.2S)|104, LOS, |33|l 




collision of two fronts (4.8) 
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